Quantum computing with magnetically interacting atoms 
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We propose a scalable quantum-computing architecture based on cold atoms confined to sites of a 
tight optical lattice. The lattice is placed in a non-uniform magnetic field and the resulting Zeeman 
sublevels define qubit states. Microwave pulses tuned to space-dependent resonant frequencies are 
used for individual addressing. The atoms interact via magnetic-dipole interactions allowing imple- 
mentation of a universal controlled-NOT gate. The resulting gate operation times for alkalis are on 
the order of milliseconds, much faster then the anticipated decoherence times. Single qubit opera- 
tions take about 10 microseconds. Analysis of motional decoherence due to NOT operations is given. 
We also comment on the improved feasibility of the proposed architecture with complex open-shell 
atoms, such as Cr, Eu and metastable alkaline-earth atoms with larger magnetic moments. 

PACS numbers: 03.67.Lx,32.80.Pj 



I. INTRODUCTION 

Over the last decade, a variety of architectures for 
quantum computing (QC) has been proposed [lj. In par- 
ticular there is a number of proposals based on neutral 
atoms and molecules trapped in optical lattices. These 
proposals focus on various realizations of the multi-qubit 
logic such as Rydberg-atom gates 0, controlled colli- 
sions y, lJJ3 electrostatic interaction of heteronuclear 
molecules 6], etc (see a popular review J2J-) While 
there is a variety of approaches available, the technolog- 
ical difficulties so far impede practical implementation 
of these schemes. Here we propose a scalable quantum- 
computing architecture which further builds upon the 
well-established techniques of atom trapping. Compared 
to thepopular neutral-atom QC scheme with Rydberg 
gates , the distinct features of the present proposal are: 
(i) individual addressing of atoms confined to sites of ID 
optical lattice with unfocused beams of microwave radia- 
tion, (ii) coherent "always-on" magnetic-dipolar interac- 
tions between the atoms, and (iii) substantial decoupling 
of the motional and inner degrees of freedom. The Hamil- 
tonian of our system is identical to that of the QC based 
on nuclear magnetic resonance (NMR) (8j, and already 
designed algorithms can be adopted for carrying compu- 
tations with our quantum processor. An implementation 
of the celebrated Shor's algorithm with a linear array of 
qubits has been discussed recently Q- 

One of the challenges in choosing the physical sys- 
tem suitable for QC is the strength of the interparti- 
cle interaction. Before proceeding with the main dis- 
cussion, let us elaborate on the suitability of magnetic- 
dipolar atom-atom interactions for QC. Compared to 
the dominant electrostatic interaction between a pair 
of atoms, magnetic-dipole interaction between a pair of 
atoms is weak (it is suppressed by a relativistic factor of 
a 2 s=s (1/137) 2 .) Yet the dipolar interaction exhibits a 
pronounced anisotropic character: the strength and the 
sign of the interaction depend on a mutual orientation 
of magnetic moments of the two atoms. Namely the 
anisotropy of the interaction plays a decisive role in realiz- 



ing a universal element of quantum logic: the two-qubit 
controlled-NOT (CNOT) gate [13 . As to the strength 
of the interaction, it determines how fast the two-qubit 
gate operates. As the interaction strength decreases, 
the operation time, tcnot, increases. Still, quantita- 
tively, the operation of the gate must be much faster 
than decoherence. For atoms in far-off-resonance optical 
lattices, the decoherence times for internal (hyperfine) 
states, chosen as qubit states, are anticipated to be in 
the order of minutes pUj. On the other hand, we show 
that for magnetically-interacting atoms placed in tight 
optical lattices tcnot is a few milli-second long. Thus, 
although the interaction is weak, it is still strong enough 
to allow for more than 10 operations on a pair of qubits. 
Considering that trapping millions of atoms is common 
now, the scalable quantum computer (QC) proposed here 
may present a competitive alternative to other architec- 
tures. 

According to DiVincenzo 12], the physical implemen- 
tation of QC should satisfy the following five criteria: 
(i) A scalable physical system with a well-characterized 
qubit; (ii) The ability to initialize the state of the qubit 
to a simple fiducial state; (iii) A "universal" set of quan- 
tum gates; (iv) Long relevant decoherence times, much 
longer than the gate operation time; (v) A qubit-specific 
measurement capability. Below, while describing our pro- 
posed QC, we explicitly address these criteria. We also 
compare it with QC based on nuclear magnetic resonance 
(NMR) techniques Q , because of the equivalence of mul- 
tiparticle Hamiltonians of our QC and the NMR systems. 
Unless noted otherwise, atomic units |e| = m e = h = 1 
are used throughout; in these units the Bohr magneton 
Ms = 1/2- 



II. CHOICE OF QUBIT 

We focus on alkali-metal atoms in their respective 
ground S1/2 state, but later (Section IVIjI also comment 
on the improved feasibility of our QC for other open- 
shell ground state and metastable atoms. The atoms are 



placed in a magnetic field and the resulting Zeeman com- 
ponents define qubit states. As shown in Fig. Q] for 23 Na, 
the Zeeman energies behave non-linearly as a function 
of the magnetic field strength B. The non-linearity is 
due to an interplay between couplings of atomic elec- 
trons to the nuclear and the externally-applied fields. 
The discussion presented below can be extended to an 
arbitrary field strength, but for illustration we consider 
the high-field limit, [i B B > AE HFS /(2I + 1). Here 
A^hfs is the hyperfine splitting that ranges from 228 
MHz for 6 Li to 9193 MHz in 133 Cs and I is the nu- 
clear spin. In this regime, the proper lowest-order states 
are product states \J, Mj) \I, Mj), where Mj and Mi are 
projections of the total electronic momentum J = 1/2 
and the nuclear spin / along the B-field. We choose 
the qubit states as |1> = \J,Mj = 1/2)|I,M 7 ) and 
|0) = \J,Mj = —l/2)\I,Mi), Disregarding nuclear mo- 
ment, these states have magnetic moments of /im = /is 
and /X| ) = —fJ-B, and the associated Zeeman splitting is 
on the order of a few GHz (see Fig^). This choice allows 
driving single-qubit unitary operations with microwave 
(MW) pulses. The MW radiation has to be resonant 
with the Zeeman frequency (a*|i) — A*|o)) B. 
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FIG. 1: Zeeman effect for the ground state of 23 Na atom. The 
states are labeled with the electronic and nuclear magnetic 
quantum numbers Mj and Mi in the strong- field limit. A 
possible choice of qubit states is shown. 

As to the initialization of the qubits of the proposed 
QC, the conventional techniques of optical pumping may 
be used. State-selective ionization may be employed for 
the read-out of the results of calculations. Ion optics 
for registering the final states is described in Ref. [g. 
Both the initialization process and the projective readout 
favor the proposed QC in comparison to the liquid state 
NMR QC, where the ensemble averaging is required for 
read-out and relaxation of the sample is important for 
initialization. 



III. INDIVIDUAL ADDRESSABILITY 

In NMR [3, the nuclear spins (qubits) are distin- 
guished by their different chemical shifts affecting reso- 



nance spin- flip frequencies. Here to individually address 
the atoms, we confine ultracold atoms to sites of a one- 
dimensional optical lattice and introduce a gradient of 
magnetic field, so that the Zeeman frequency depends on 
the position of the atom in the lattice (see Fig|5J). This 
addressing scheme was investigated in details in Ref. [13j , 
in the context of quantum processor using trapped ions. 
It is also worth mentioning a similar idea for QC with 
heteronuclear diatomic molecules 6] , where a gradient of 
electric field was applied along optical lattice affecting 
flip frequencies of molecular dipoles. 




FIG. 2: Proposed quantum computing architecture. The 
atoms are confined to sites of an optical lattice. Non-uniform 
magnetic field is required for individual addressing of the 
atoms with pulses of the resonant microwave radiation. 



A. Optical lattice 

The ID potential of the optical lattice created by 
a standing wave CW laser beam reads V op t(z) — 
Vbsin 2 (^fz). The depth of the wells is Vq = 8na x 
a a (u>)Ii,, where a n (wi) is the dynamic electric-dipole po- 
larizability of the atom, II is the laser intensity, and 
a « 1/137 is the fine-structure constant. Depending on 
the detuning of the laser frequency u>l from a position 
of atomic resonance, the polarizability can accept a wide 
range of values. To restrict the transverse atomic motion 
in a ID optical lattice one must require that a a (uJL) > 0. 
The case of negative a a (u>L), although requiring 3D op- 
tical lattice, offers an advantage of reduced photon scat- 
tering rates. Since the atoms would be located at the in- 
tensity minima, the rates (in a tight confinement regime) 
would be suppressed by a factor of 1/2^/Er/Vo <C 1, 
where En = a 2 w|/(2Af) is the photon recoil energy for 
an atom of mass M. Further, an atom is assumed to oc- 
cupy the ground motional state of the lattice wells and 
one atom per site filling ratio is assumed. Techniques for 
loading optical lattices are being perfected [lj, [l5j . 

In the ID optical lattice, the neighboring atoms are 
separated by a distance of R = A/2. To maximize the 
magnetic-dipole interaction (ex 1/-R 3 ) between two neigh- 
bors we require that the frequency of the laser lul is cho- 
sen as high as possible. A natural limit on lu^ is the 
ionization potential (IP) which ranges from 3.8 eV for 
Cs to 5.4 eV for Li. Practically, ujl must be somewhat 



below the IP to avoid resonances in the high-density of 
states near the continuum limit. In the estimates below 
we use ujl ~ 5 eV. This corresponds to 



A, 



250 nm = 4700 bohr . 



values of magnetic moments for the two qubit states (see 
Section lV Aj l. For the parameters chosen above the mag- 
netic forces are several orders of magnitude smaller than 
the optical ones. 



While working with such short ultra-violet wavelengths 
seems feasible, more conventional wavelenghts of 400 — 
600nm are adequate when working with complex open- 
shell atoms with larger magnetic moments (see Sec- 
tion ED. 

One additional requirement is that the tunneling time 
between adjacent sites of an optical lattice is sufficiently 
long, so that the atoms maintain their distinguishability 
during the computation. The site-hoping frequency can 
be estimated as [lja 



,/ = E, 






3/4 



exp 



-2y/V /E R 



For 23 Na, the requirement that the characteristic tun- 
neling time it I J rs 10 sec leads to a barrier hight of 
Vq = 0.4 mK. At a fixed u>l, the tunneling may be 
suppressed by increasing the barrier height Vq and by 
using heavier atoms. Additionally, due to the Pauli ex- 
clusion principle, the tunneling may be suppressed by 
using fermionic atoms. 



IV. MULTI-QUBIT OPERATIONS 

Having discussed one-qubit operations, we now con- 
sider a realization of the universal controlled-NOT 
gate |lCj based on magnetic-dipolar interaction of two 
atoms. It is worth mentioning that we discuss this gate 
here only for illustrative purposes, to estimate the char- 
acteristic duration of the CNOT operation, tcnot- The 
many-body dynamics of the system with the interpar- 
ticle interaction which is "always on" is such that the 
two-qubit gates are executed as a part of the natural dy- 
namics of the system, and a s pec ial care has to taken to 
govern this natural evolution [17fl. We will elaborate on 
this point at the end of this Section. 

The basic requirement for the controlled-NOT gate [l(j 
is that the frequency of |1) — > |0) transition of the target 
atom depends on the state of the control atom. Since the 
quantization (B-field) axis is directed along the internu- 
clear separation between the two atoms, the magnctic- 
dipole interaction can be represented as 



B. B-field gradients 

The magnetic field leads to the field-dependent Zeeman 
effect and the gradient dB/dz of the B-field allows one 
to resolve the resonant frequencies of individual atoms, 
see Fig. [21 In the Section IIVI we show that a typical 
two-qubit operation has a duration tcnot of a few mil- 
liseconds. To find the gradient dB/dz we require that 
the single-qubit NOT operation, performed by the reso- 
nant 7r-pulse of microwave radiation, takes a comparable 
time tnot of 1 ms. Such pulse may be resolved by two 
neighboring atoms, provided that their resonant frequen- 
cies differ by A/ NOT = 1/(27tt NO t) «2x 10 2 Hz. The 
required field gradient is relatively weak 

dB/dz « 20 G/cm 

and is comparable to typical gradients in conventional 
magnetic traps. Much steeper gradients of 3 x 10 3 G/cm 
over a region of a few millimeters have been demonstrated 
by Vuletic et al. [Ig. These authors employed ferromag- 
netic needles that collect and focus B-fields of electro- 
magnetic coils. With such gradients the performance of 
the single- qubit gates can be substantially improved, 
t NO t ~ 10 /is. 

Another limitation on the gradients is that the mag- 
netic force fidB/dz must be much smaller than the opti- 
cal force — dV op t(z)/dz. This limitation affects not only 
the confinement but also a degree of coupling of inner and 
motional degrees of freedom because of the difference in 



V = ~W £(l + *A,0)MA(l)/i-A(2), 



(1) 



A=-l 



where [J,\(k) represents the spherical components of the 
magnetic moment operator fi for atom k. For perform- 
ing the CNOT gate operation we need to resolve the 

frequency difference of 2a 2 R~ 3 (p\i) — M|0>) ■ F° r our 
choice of parameters we find that A/cnot ~ 40 Hz. 
Although this number may seem small, it is compara- 
ble to typical coupling strengths of 20-200 Hz in NMR. 
Furthermore, the minimum duration of the MW pulse 
is tcnot = 1/(27tA/cnot) ~ 3ms, allowing for more 
than 10 4 CNOT operations during anticipated decoher- 
ence time (in the order of minutes [llj ) . 

As in the NMR, the interaction between the atoms is 
always on and special care has to be taken to stop the 
undesired time evolution of the system and carry out 
controlled calculations. Fortunately, it is straightforward 
to show that the many-particle Hamiltonian of our sys- 
tem is equivalent to that of a system of interacting spins 
in NMR and thus already developed techniques can be 
adopted. Specifically, the Hamiltonian is 



-HnMR = Y^ Ui ( a ^i + o Y 9l 3 (°*)* (°*) 



2 *—>■ 



(2) 



Here (p^) i are the Pauli matrices for an atom i located 
at position Zj. Introducing the average magnetic moment 
p, = (yUm+/U|o))/2 and the difference S /i = (/xiu — jj,\q\)/2, 



the one-particle frequencies may be expressed as 



u>i = -6niB(zi) + ^ ' 

3 

and the two-body couplings as 



2a 2 



TfijSlM 



gij = -2a 2 /(\z t - Zj\ 3 ) S^d/Xj . 

Since the Hamiltonian Hnmr, Eq.Q, is identical to 
that of the NMR based QC [l^j, the already developed 
NMR algorithms can be adopted. For example, the time 
evolution for a given atom may be reversed by invert- 
ing populations. The same idea is applicable for the in- 
terparticle couplings. The CNOT gate may be imple- 
mented by applying a sequence of one-qubit transforma- 
tions (short pulses) and allowing a given coupling to de- 
velop in time | l7| . An implementation of the celebrated 
Shor's algorithm on the system of linear array of qubits 
has been discussed recently by Fowler et al. 0]. 

An important point is that the two-particle couplings 
result in a coherent development of the system. As we 
demonstrate below, the short NOT pulses introduce neg- 
ligible decoherences due to excitation of motional quanta 
during the gate operation. 



V. DECOHERENCES 

Finally, we address possible sources of decoherences. 
We divide these sources into two classes: (i) induced due 
to the required operation of the quantum processor (due 
to performing gates) and (ii) architectural, i.e., due to 
traditional sources of decoherences, such as photon scat- 
tering (e.g. Raman), instabilities of laser output, insta- 
bilities in the magnetic field, collisions with background 
gas, etc. Below we analyze the decoherences caused by 
the gate operations and we also estimate the upper limit 
on tolerable fluctuations of the magnetic field. The re- 
maining decoherences depend on the details of experi- 
mental design, e.g., laser wavelengths available, and we 
can not fully address them at this point. 



A. Motional heating 

The atomic center-of-mass (CM.) evolves in a com- 
bined potential of the optical lattice and magnetic field. 
This potential depends on the internal state of the atom, 
due to a difference in dynamic polarizabilities and mag- 
netic moments of the states. When population is trans- 
ferred from one qubit (internal) state to the other, the 
motion of the CM. is perturbed and the atom may leave 
the ground state of CM. potential (motional heating). 
Here we calculate the relevant probability and show that 
it is negligible. We find that the motional heating is 
suppressed due to adiabaticity of the population-transfer 
process. In other words, the inner and motional degrees 



of freedom decouple because our quantum processor is 
relatively slow. It is worth emphasizing that in the pop- 
ular neutral-atom QC proposal |2| , the operation of the 
gates is much faster and the motional heating is of con- 
cern [is|. 

As discussed in the previous section, the QC proposed 
here requires only single-qubit rotation pulses with a 
characteristic duration tnot- The two-qubit gate op- 
erations are performed by the NOT pulses and due to 
the natural coherent dynamics of the system, so it suffi- 
cient to consider the coupling of the inner and motional 
degrees of freedom due to the NOT pulse only. 

For a resonant radiation, during the MW pulse, the 
effective Hamiltonian may be represented as 

H(R,r) = ^(|0)(1| + |1)(0|)+ (3) 

^ + %(R)|1)(1| + V| 0) (R)|0)(0|. 

Here R and P are the CM. coordinate and momentum, 
and r encapsulates internal degrees of freedom (|0) and 
|1)). The internal states are coupled via Q = vt/tnoT) 
the Rabi frequency of the MW transition. The CM. 
evolves either in Vj ) or Vm potentials. We treat the 
difference between the CM. potentials W — Vm — Vj \ 
as a perturbation and characterize the CM. motion with 
eigenstates |$ n ) and energies S n computed in the Vj ) 
potential. 

/ P 2 
\2M 



V|o)(R) |*»)=£»|*r. 



We assume that at t = the atom is in the motional 
and internal ground state ^(R, r,< = 0) = |0) |3>o)< As 
a result of the transfer of population (NOT operation, 
|0) — ► |1) ) there is a finite probability P p for an atom to 
end up in one of the excited motional eigenstates |$ p ). 

We evaluate the probability P p perturbatively by ex- 
panding 

*(R,r,t) « o(«)|0) |$o> +b(t)\l) |$ > + c (t)|l) |$ p ). 

Substituting this expansion into the Schrodinger equa- 
tion and solving it perturbatively (perturbation is W, the 
two-level Rabi evolution is treated exactly ) we obtain 



a(t) 

b(t) 



. d , , 



sin(Q/2i) , 

cos{tt/2t) , 

sin (n/2t) exp (iuj p0 t) (<5> P \W\<£> ) 



where uj p o = £ p — £o- Calculation of the probability based 
on c(i) requires certain care, since after the initial appli- 
cation of the pulse the perturbation W remains on in- 
definitely Using an approach discussed in Ref. [lflj . we 
arrive at the probability of excitation 



-P p (tnot) 



(*p|W|4o) 



w p o 



e(2a. 



(4) 



where we introduced the adiabaticity function 



G(0 



e 



ie 



i) 



l + £ 2 -2£sin( 



2f 



(5) 



This function is plotted in Fig. [3J The function is 
bounded G (£) < 1; there is no singularity at the res- 
onant frequency G(£ = 1) sw 0.87, since the CM. experi- 
ences only a quarter of the period of oscillating function 
during the 7r pulse. Let us investigate various limits of 
G (£) . For an instantaneous turn-on of the perturbation, 
G(£ »1)«1, so that 



P P (tnot) 



(*p|W|*o) 



UJpO 



TNOT^pQ <C 1 . 



For a slow (adiabatic) application of the 7r-pulse G(£ <C 
1) rj l/£ 2 . In this adiabatic limit 



-P p (tnot) 



($ P |VF|$o) 



n/2 



, TNOT^pQ ^> 1 . 



the dynamic polarizability a Q (w) for the two qubit 
states. This d ifference leads to a modification of 
Who oc -\/ a a(w). The relevant probability may be 
estimated as ^ ~ G (fi/4who) 55 {^^bB/Alo) s=s 
10~ 10 (2/j,bB/Au>) , where Aw is a detuning of the laser 
frequency from the energy of the intermediate state con- 
tributing the most to the dynamic polarizability a a (oJL)- 
(Here we used the field gradient dB/dz — 20G/cm, so 
that tnot ~ tcnot, see section IIIIBI ) In practice, 
2fi B B/Auj < 1, so that P 2 < 10~ 10 . This probabil- 
ity can be reduced even further by adjusting the laser 
wavelength so that the dynamic polarizabilities of the 
two qubit states are the same [l8j . 

Finally, we would like to emphasize that the motional 
heating is not an issue for our processor, because our QC 
is relatively slow. The CM. cycles through many oscil- 
lations while the perturbing MW field is applied. This 
adiabatic averaging leads to the suppression of the mo- 
tional heating. 




FIG. 3: Adiabaticity function G(£), Eq.©. 

Now we can carry out the numerical estimates. To 
compute the matrix elements of the perturbation, we ap- 
proximate the bottom of the optical potential as that of 
harmonic oscillator of the frequency Who = t 2 - \/2Vq/M 
and use the harmonic oscillator wavefunctions for |$ n ). 
For 23 Na and our lattice parameters (Aj, « 250 nm, Vq rs 
0.4 mK ), w^ ~ 2?r x 2 x 10 6 Hz, while the Rabi fre- 
quency (for the extreme case of dB/dz — 2 x 10 3 G/cm) 
Q/2 = 2tt x 0.25 x 10 5 Hz, i.e. even in the case of 
steep gradients we deal with an adiabatically slow pulse 
and the motional heating is suppressed by a factor of 
« [w p0 /(tt/2)] 2 . 

For a perturbation due to an interaction of mag- 
netic moments with the gradient of magnetic field W = 
(/i|o) — Mil)) {dB/dz)z. For our parameters, the proba- 
bility of exciting motional quantum is just P\ ~ 6 x 10~ 7 . 

Another source of coupling of internal and mo- 
tional degrees of freedom is due to the difference in 



B. Decoherences due to fluctuations of the 
magnetic field 

The atoms in our quantum processor are required to 
have magnetic moments. This magnetic moments can 
couple to the ambient magnetic field. The fluctuations 
of these ambient fields can cause the qubit states to lose 
coherence. In addition, the Johnson/shot noise in the 
coils providing the gradient of the B-field can lead to the 
decoherences as well. Here we estimate the upper limit 
of tolerance to these fluctuations. 

Before we carry out the estimate, we notice that the 
sensitivity to fluctuations of the magnetic fields can be 
avoided altogether using so-called Decoherence Free Sub- 
spaces (DFS) mmEa. The idea of the DFS is to rede- 
fine the qubit states using linear combinations of product 
of states of original qubits; each resulting state accumu- 
lates the same phase due to environmental interactions. 
Since the total phase of wavefunction is irrelevant, the 
DFS is highly stable with the respect to the external 
perturbations. This powerful DFS approach has been 
already experimentally verified in a number of cases, in- 
cluding trapped ions |23|. Nevertheless, an introduction 
of the DFS requires reconsideration of the NMR algo- 
rithms, which is beyond the scope of the present paper. 
Below we estimate the decoherence of our original qubit 
defined as a couple of magnetic substates of an atom. 

To make an order-of-magnitude estimate, we assume 
that the fluctuating field B'(t) is along the quantization 
axis and it has a white noise spectrum. With the fluctu- 
ating field present, the accumulated phase difference be- 
tween the qubit states is 8(f) ~ 2 f Q h\q)B' (ti) dt\. (Here 
we assumed that /Urn ~ — /^lo) ■) Averaging over different 



realizations 



(expW)) « 1 + i <J0) - - <^ 2 > 



1 + 2*/ a*|o> (B' (ti)> dt 



-2/i : 



o) 



o Jo 



B' (h) B' (t 2 )) dhdh . 



This expression may be simplified using ensemble aver- 
age (B' (ti)) = 0, and the autocorrelation function for a 
stochastic process with no memory 

(B'(t 1 )B'(t 2 )) = (B 2 )j(ti-t 2 ), 

where (-B 2 ) = const. If a measurement is carried out 
after tcnot, the probabilities would differ from the exact 
result by 

e B = ((cxp (i54>)) - l) 2 sa 4/i 2 0> (B 2 )^ t cnot ■ 

According to Knill |2J|, this error can be as high as 1%. 
This leads to a tolerable level of the B-field noise 



T 



(£ 2 L £ 10 ~ 10 



This limit is relatively easy to attain experimentally |25| . 

To summarize the main results of this section, we have 
demonstrated that for our quantum processor the mo- 
tional heating caused by the gate operations is negligible. 
The derived formula Q is also applicable for an analysis 
of coupling of inner and motional degrees of freedom of 
QC based on heteronuclear molecules p. We also have 
evaluated the tolerable level of fluctuations in magnetic 
field. It is worth emphasizing that the sensitivity to en- 
vironmental field noise can be greatly eliminated using 
decoherence free subspaces, as discussed in the beginning 
of this subsection. 

There are other potential sources of decoherences, e.g. 
photon scattering, instabilities of laser output, and colli- 
sions with background gas. However, we do not address 
them at this point since they depend on specifics of ex- 
perimental design. Still it is anticipated |ll| that the 
decoherence times for internal states may be in the order 
of minutes. 



VI. COMPLEX OPEN-SHELL ATOMS 



For example, recently cooled and trapped Cr and Eu 
atoms 26j have magnetic moments of stretched states 
MCr = 6/ib and fiEu — 7/ib- Thus if a laser of a 
400 nm wavelength is employed for constructing the op- 
tical lattice, the CNOT-gate operations would require 
tcnot ~ 0.1ms. The use of 30-mW violet 400 nm LD 
laser has been recently demonstrated in cooling and trap- 
ping experiments |27| . The requirement tcnot ~ 1 ms 
would bring the wavelength to an even more practical 
800 nm. 

Another alternative is to use long-lived [28| metastable 
divalent atoms, Mg, Ca, Sr and Yb. Several groups have 
demonstrated cooling and trapping these atoms 29] , with 
Yb BEC recently attained 30] . Here the state of interest 
is 3 P2, with [j,3p 2 = 3/is and for a 400 nm laser tcnot ~ 
lms. 



VII. CONCLUSION 

In summary, we proposed a scalable quantum comput- 
ing architecture based on magnetically-interacting cold 
atoms confined to sites of a tight optical lattice. The 
atoms are placed in a non-uniform magnetic field and the 
individual addressing is attained by pulses of microwave 
radiation for a minimum duration of 10 /is. The uni- 
versal two-qubit CNOT gates require times in the order 
of milliseconds. The multi-particle Hamiltonian of the 
system is equivalent to that of the QC based on NMR 
techniques so the already developed quantum algorithms 
may be adopted. 

Compared to the popular neutral-atom QC scheme J2J, 
the distinct advantages of the present proposal are: (i) in- 
dividual addressability of atoms with unfocused beams of 
microwave radiation, and (ii) coherent evolution due to 
the "always-on" magnetic-dipolar interactions between 
the atoms, (iii) substantial decoupling of the motional 
and inner degrees of freedom. While the main disadvan- 
tage of our QC is the slow rate of operation, it is antici- 
pated that one could carry out 10 CNOT operations on 
a single pair of atoms before the coherence is lost. When 
addressing tens of thousands of atoms is done in paral- 
lel, this would amass 10 8 two-qubit operations and 10 10 
single-qubit operations. 



In this paper we focused on alkali-metal atoms (/i = 
Hb) and we found that wavelengths of A ~ 250 nm are 
needed for tcnot ~ 3 ms. While working with such short 
ultra-violet wavelengths still seems feasible, the restric- 
tion may be relaxed by using complex open-shell atoms 
with larger values of magnetic moments. This allows in- 
creasing the interaction strength (ex /z ) and thus working 
with more conventional laser wavelengths. 
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